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NOTATION 

Aft  perpendicular 
Buoyancy 

Residual  resistance  coefficient 
Frictional  resistance  coefficient 
Maximum  diameter  of  body 
Stability  index 

Moment  of  inertia  about  y and  z axes,  respectively 
Length  of  body 

Distance  from  nose  to  center  of  buoyancy 
Mass  of  body 

Moment  about  y and  z axes, respectively 
Parallel  middlebody 
Perturbation  angular  velocity  about  y and  z axes,  respectively 
Total  surface  area  of  body 
Time 

Perturbation  velocities  in  the  x,y,  and  z directions, 
respectively 

Equilibrium  forward  speed 
Forward  speed  in  knots 
Weight 

Right-handed  Cartesian  coordinates  with  origin  fixed  on 
the  body  axis  at  the  center  of  buoyancy;  positive  x is 
directed  forward  along  the  longitudinal  body  axis 
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XB*  XG 


ZB’  ZG 


Longitudinal  coordinates  of  centers  of  buoyancy  and 
gravity,  respectively,  relative  to  x,y,z  coordinate  system 

Force  in  the  y direction 

Vertical  coordinates  of  centers  of  buoyancy  and  gravity , 
respectively,  relative  to  x,y,z  coordinate  system. 


Force  in  the  z direction 

Perturbation  angular  displacement  about  y axis 
Equilibrium  pitch  angle  relative  to  the  horizontal 
Kinematic  viscosity  of  water 
Mass  density  of  water 
Nondimensional  stability  root 
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The  stern  control  surfaces  of  a 60-foot  version  of  the 
DOLLY  VARDEN  Buoyant  Rise  Test  Vehicle  were  designed  from  the 
standpoint  of  dynamic  stability  in  vertical  flight.  Computer 
predictions  were  made  of  the  vertical  flight  trajectories  of  the 
vehicle.  It  was  found  that  a cruciform  configuration  of 
all-movable,  rectangular  stern  appendages  with  tip  to  tip  span 
of  10  feet  and  chord  length  of  5 feet  will  provide  adequate 
stability  in  vertical  buoyant  rise  flight  and  required  maneuverability 
for  recovery.  It  was  also  determined  that  a net  buoyancy  of 
20000  pounds  would  provide  the  required  operational  speed  in 
vertical  buoyant  rise  flight. 


ADMINISTRATIVE  INFORMATION 

Predictions  of  the  motions  of  the  DOLLY  VARDEN  Buoyant  Rise  Test 
Vehicle  were  made  for  the  Ship  Acoustics  Department  of  the  David  W.  Taylor 
Naval  Ship  R&D  Center  in  accordance  with  memorandum  1946:JS:DN  of 
13  February  1976. 


INTRODUCTION 


The  Ship  Performance  Department  of  the  David  W.  Taylor  Naval  Ship  R&D 
Center  (DTNSRDC)  was  requested  by  the  Ship  Acoustics  Department  to  provide 
assistance  in  developing  the  general  arrangements,  designing  the  control 
surfaces,  and  predicting  the  motions  and  trajectories  of  the  DOLLY  VARDEN, 
Buoyant  Rise  Vehicle,  a newly  designed  acoustic  test  vehicle  similar  to 
the  KAMLOOPS.  Although  the  DOLLY  VARDEN  may  eventually  operate  in  powered, 
submerged  horizontal  flight,  it  will  be  initially  used  as  a buoyant  rise 
vehicle,  and  will  be  required  to  ascend  with  the  longitudinal  axis  in  a 
nearly  vertical  orientation  or  attitude. 

It  is  intended  that  the  vehicle  be  given  a specified  net  buoyancy  and 
released  from  far  below  the  surface  with  an  initially  vertical  attitude.  The 
vehicle  will  accelerate  upward  rapidly  and  must  attain  a specified  speed  of 
40  knots  at  a depth  of  400  feet.  At  200  feet  below  the  surface,  a recovery 
maneuver  is  to  be  initiated.  It  is  anticipated  that  the  stemplanes  would 
be  deflected,  thus  turning  the  vehicle  from  a nearly  vertical  flight  to  a 
more  nearly  horizontal  flight  path  with  a consequent  slowing  of  its  forward 
speed.  During  the  vertical  portion  of  the  flight,  it  is  also  required  that 
deviation  in  the  vehicle  attitude  from  the  vertical  and  the  horizontal 
displacement  from  the  launch  point  not  be  excessive. 

An  analysis  of  the  dynamic  behavior  of  the  DOLLY  VARDEN  in  vertical 
buoyant  rise  flight  is  presented  in  this  report.  Hydrodynamic  coefficients 
for  the  buoyant  rise  vehicle  have  been  estimated,  and  the  vertical  flight 
dynamic  stability  has  been  examined  as  a function  of  net  buoyancy  and 
center  of  gravity  location.  A set  of  stem  appendages  is  recommended  to 
provide  acceptable  stability  and  maneuverability. 
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The  motion,  trajectory,  and  speed  of  the  vehicle  during  the  vertical 
flight  and  proposed  recovery  maneuver  have  been  determined  using  the  sub- 
marine motion  computer  simulation  program  ZZMN.  The  effect  of  amount  of 
net  buoyancy  is  examined.  In  addition,  the  influence  of  the  maximum  stern- 
plane  angle  and  rate  of  deflection  on  the  recovery  maneuver  is  investigated. 

GEOMETRIC  CHARACTERISTICS 

The  DOLLY  VARDEN  is  a large  test  vehicle  roughly  similar  in  external 
appearance  to  the  Navy's  modern  attack  submarines.  The  vehicle  is  a body 
of  revolution  with  typical  forebody,  parallel  middlebody,  and  afterbody, 
and  with  cruciform  tail  appendages  the  size  and  shape  of  which  are 
dependent  upon  stability  and  maneuverability  considerations.  For  the 
buoyant  rise  configuration,  there  is  no  propeller  or  bridge  fairwater. 

The  length  of  the  vehicle  is  approximately  60  feet  with  a 23-foot 
parallel  middlebody  section.  In  the  future  the  parallel  middlebody  may 
be  extended  to  accommodate  a propulsion  unit.  However,  it  should  be 
emphasized  that  the  analysis  presented  here  will  apply  only  to  the  60-foot 
version. 

Figure  1 shews  the  outline  of  the  bare  hull  of  DOLLY  VARDEN.  Dimen- 
sions are  as  provided  by  the  DTNSRDC  Ship  Acoustics  Department.  Table  1 
shows  various  geometrical  characteristics  of  the  bare  hull. 

The  weight  and  longitudinal  location  of  the  center  of  gravity  are 
treated  as  independent  variables  in  the  analysis.  In  the  vertical  attitude, 
the  desired  upward  acceleration  of  the  vehicle  is  produced  by  a net  dif- 
ference between  the  buoyancy  and  weight,  the  magnitude  of  which  affects 
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Table  1 - Geometric  Characteristics  of  the 
DOLLY  VARDEN  (Bare  Hull) 


Characteristic 

Value 

Length 

59.42  feet 

Diameter 

(Maximum) 

7 . 6B57  feet 

Length:  Nose 

PMB 

Afterbody 

14 .090  feet 
22 .972  feet 
22.358  feet 

Volume 

2125  feet3 

Wetted  Surface 
Area 

1195  feet2 

Longitud ina  1 
Center  of  Buoyancy 
(from  nose) 

27.9  feet 

Vertical  Center 
of  Buoyancy 

0.0  feet 

Mass  , 

Neutral  Buoyancy 

4.120  X 103  ^ 
s lugs 

Buoyancy 

1 .326  X 105 
pounds 

Longitudinal 
Moment  of  Inertia 

0.75  X 106 

slug-feet 

* Mass  density  of  freshwater  at  4 degrees  C.  (^  = 1.937  slugs/feet3) 
**  Acceleration  due  to  gravity  (g  =.  32.175  feet/  sec2) 


both  speed  and  stability.  The  center  of  gravity  is  generally  located 
aft  of  the  center  of  buoyancy  in  order  to  provide  stability  in  vertical 
flight.  In  the  analysis,  the  net  upward  buoyancy  ranges  from  16,000 
to  24,000  pounds  and  the  center  of  gravity  is  located  1.5  to  2.5  feet  aft 
of  the  center  of  buoyancy. 

The  stability  characteristics  of  a neutrally  buoyant  configuration 
with  centers  of  buoyancy  and  gravity  coincident  are  also  investigated, 
since  this  condition  might  be  obtained  for  a zero  trim,  neutrally  buoyant 
vehicle  in  powered  horizontal  flight. 

The  center  of  gravity  of  a submarine  would  normally  be  displaced  from 
the  axis  of  symmetry  in  order  to  provide  static  stability  in  roll  and  pitch. 
In  the  present  analysis,  however,  the  effect  of  this  displacement  of  the 
center  of  gravity  has  not  been  considered  and  it  is  assumed  to  lie  on  the 
body  longitudinal  axis  of  symmetry  in  all  cases. 

The  moment  of  inertia  of  the  vehicle  in  pitch  depends  on  the  ballast 
conditions.  However,  the  differences  in  the  actual  ballast  conditions 
investigated  occur  near  the  middle  of  the  vehicle,  and  thus  their  effect 
on  the  moment  of  inertia  is  minimal.  The  value  of  the  moment  of  inertia 
used  in  the  analysis  was  computed  for  the  16 ,000  pound  net  buoyant  vehicle 
with  uniform  mass  distribution,  and  was  assumed  to  be  the  came  for  all 
other  ballast  conditions. 

DYNAMIC  STABILITY  AND  DESIGN  APPENDAGES 

The  sizes  of  the  stern  appendages  for  the  DOLLY  VARDEN  are  based 
upon  requirements  for  adequate  stability  and  maneuverability  of  the  vehicle 
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in  vertical  buoyant  rise  flight.  The  linear,  vertical  plane  equations  of 
motion  are  given  in  Appendix  A,  and  the  dynamic  stability  characteristics 
of  the  vehicle  configurations  having  a range  of  appendage  chord  lengths  and 
ballast  conditions  are  obtained  from  the  roots  of  the  governing  character- 
istic equations.  It  is  assumed  that  the  body  is  initially  in  steady  straight- 
line  flight,  moving  in  a purely  forward  direction  with  speed  U at  an  angle 
0q  relative  to  the  horizontal  as  shown  in  Figure  2.  The  coupled  normal 
force  and  pitch  moment  equations  for  small  perturbations  from  the  equili- 
brium  condition  of  steady  flight  are  given  by  Equations  (A-l)  and  (A-2)  of 
Appendix  A,  respectively,  and  the  characteristic  equation  corresponding  to 
the  vertical  plane  motion  is  a cubic  of  the  form 

ac rj  + b <rKz  + c<rK  + d = o (D  j 

where  (T  is  the  nondimensional  root,  and  the  coefficients  are  defined  by 
Equations  (A-6)  through  (A-9) . 

Of  particular  interest  is  the  condition  where  the  vehicle  is  moving 
vertically  upward.  The  metacentric  derivatives  ' and  Mg'  given  by 
Equations  (A-3)  and  (A-4) , respectively,  for  0^  *>  90  degrees  are  evaluated 
in  Table  2 for  several  combinations  of  the  net  buoyancy  and  longitudinal 
center  of  gravity,  x_.  Since  the  nondimensional  metacentric  derivatives 

(j 

are  speed  dependent,  the  values  indicated  in  Table  2 have  been  multiplied 
by  the  square  of  the  speed  in  knots.  A negative  x„  indicates  that  the 
center  of  gravity  is  aft  of  the  center  of  buoyancy. 


Fipure  2 


- Definition  Sketch  of  Problem  Geometry 
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Table  2 - Ballast  Conditions  and  Metacentric 
Derivatives 


Net 

Bnovancy 
in  pounds 

We ight 
in  pounds i 

Mass 

z«vk  * 

in  knots'1 

XG 

in  feet 

— 

* 

in  knotsz 

16000 

1 .166  X 105 

0.01782 

1 .6385 

-1  .5 
-2.0 
-2.5 

-0 .3014 
-0.4019 
-0.5023 

20000 

1.126  X 105 

0.01721 

2.0482 

-1 .5 
-2  .0 
-2.5 

-0.2911 
-0  .3881 
-0.4852 

24000 

1 .086  X 105 

0.01659 

2.4578 

-1  .5 
-2.0 
-2.5 

-0.2807 

-0.3743 

-0.4678 

J 

_ l 

* 


Q0-  90  Degrees 


r 
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A cruciform  appendage  conf iguration  with  tip  to  tip  span  of  10  feet 
was  specified,  and  for  simplicity  of  construction,  all-movable,  rectangular 
appendages  were  preferred.  The  trailing  edges  of  the  appendages  were  fixed 
at  56.6  feet  aft  of  the  nose  and  the  effect  of  appendage  chord  length 
variation  on  the  vehicle  dynamic  stability  was  investigated  over  the  range 
of  ballast  conditions  indicated  in  Table  2. 

The  hydrodynamic  stability  derivatives  for  vehicle  configurations 
with  3,  4,  5,  6,  7 and  8-foot  chord  appendages  have  been  estimated.  These 
stability  derivatives  in  nondimens ional  form  are  given  in  Table  3 and 
follow  the  notations  and  conventions  of  Reference  1.  The  DOLLY  VARDEN 
is  geometrically  symmetric  relative  to  the  longitudinal  axis,  hence  the 
values  for  the  vertical  plane  derivatives  also  correspond  to  the  hori- 
zontal plane  derivatives  indicated  in  parentheses.  A sketch  of  the  5-foot 
chord  appendage  is  shown  in  Figure  3. 

i 

The  roots  of  the  characteristic  equation  <T^  determine  the  nature  of 

the  vertical  plane  motion.  For  the  vehicle  configuration  with  5-foot 

chord  appendages,  the  characteristic  equations  corresponding  to  the  ballast 

conditions  of  Table  2 are  provided  in  Table  4,  where  V is  the  speed 

K. 

in  knots  and  the  appropriate  metacentric  and  hydrodynamic  derivatives 

are  given  from  Tables  2 and  3,  respectively.  The  equations  in  Table  4 

have  been  divided  by  the  coefficient  of  the  cubic  term  A which  is 

positive  and  relatively  unchanged  in  all  cases. 

/ 

The  nondimens ional  stability  roots  Cfc  are  shown  in  Figure  4 for 

the  5-foot  chord  configuration  with  net  buoyancy  of  16,000  pounds,  and 

/ 

center  of  gravity  = -1*5  feet.  The  symbol  Re  denotes  the  real  roots. 
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Table  3 - Stability  Derivatives  for  the 
DOLLY  VARDEN 


Chord  Length  in  feet 


Coefficient ' 

3 

4 

in 

6 

7 

8 

zw 

(y'v) 

-.02575 

- .02814 

-.02973 

-.03083 

- .03159 

-.03215 

Mw 

(-< ) 

+.006814 

+.006096 

+.005664 

+.005483 

+.005466 

+ .0055r,3 

(-Yf) 

- .01111 

-.01182 

-.01226 

-.01244 

-.01246 

-.01237 

Mq 

(n') 

-.005267 

-.005454 

-.005516 

-.005477 

-.005369 

-.005217 

Z* 

<▼*) 

- .02068 



> 

-9 

f 

K 

(-<) 

0.0 

> 

> 

> 

> 

7\ 

(-Y*  ) 

0.0 

> 

, 

— - 

V 

(Nf) 

- .001473 

- 

> 

> 

** 

> 

* 

For  definition  of  nondimensional  coefficient,  see  Reference  1. 
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Longitudinal  Axis  in  feet 

Figure  3 - Proposed  Stern  Appendage  for  the  DOLLY  VARDEN 


Vertical  Stability  Roots 
Rea  1 Root  (Re 


Imaginary  Root  (Imd*) 
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Re  (f* , denotes  the  real  part  of  the  complex  root,  and  Im  g"  , denotes 
the  Imaginary  part  of  the  complex  root.  For  speeds  less  than  15  knots, 
there  Is  a positive  real  root  indicating  a divergent  path  and  a pair  of 
complex  roots  with  negative  real  part  indicating  that  the  response  is 
oscillatory.  For  speeds  above  15  knots,  all  of  the  roots  are  real,  and 
one  of  the  roots  is  positive  indicating  a nonoscillatory  divergent  response 
for  the  vehicle.  The  small  positive  real  root  <j^  means  that  the  pitch 
motion  is  unstable  at  all  speeds  for  the  particular  appendage  and  ballast 
condition  indicated.  At  higher  speeds,  the  positive  root  approaches  zero. 

It  will  be  seen  later  that  the  small  positive  root  plays  a minor  role  in 
the  trajectories  of  the  vehicle.  The  stability  analysis  is  not  strictly 
valid  in  the  very  low  speed  range  where  the  vehicle  is  accelerating 
rapidly.  Simulation  of  the  trajectories  indicates  the  vehicle  is  actually 
stable  at  and  immediately  following  launch,  but  then  has  essentially  the 
stability  characteristics  shown  in  Figure  A as  the  vertical  acceleration 
decreases . 

The  DOLLY  VARDEN  will  be  operating  primarily  in  the  high  speed  range, 

/ 

and  the  very  small  value  of  the  unstable  root  Re  o'  is  not  expected 
to  significantly  influence  the  vehicle  trajectory.  Hence,  the  effect  of 
variation  of  appendage  chord  length  and  ballast  condition  on  the  stability 
characteristics  at  infinite  speed  will  be  investigated. 

At  infinite  speed  Equation  (1)  reduces  to  a quadratic  of  the  form 

A<fKz+  B(j;+  CQ  « 0 (2) 

where  the  coefficients  are  given  by  Equations  (A-6) , (A-7),  and  (A-12)  . A 
necessary  and  sufficient  condition  for  stability  is  that  all  of  the  coefficients 
of  Equation  (3)  be  positive.  For  the  range  of  ballast  conditions  and 
chord  lengths  indicated  in  Tables  2 and  3,  respectively,  the  coefficients 
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A and  B are  positive  and  remain  relatively  unchanged  (See  Table  4). 

The  infinite  speed  stability  condition  then  requires  that  the  coefficient 
CQ  (Equation  A-12)  be  greater  than  zero.  The  infinite  speed  stability 
index  G n is  given  by  Equation  (A-14) , where  > 0 denotes  stability  and 
a larger  value  of  G^  denotes  a greater  degree  of  stability. 

In  Figure  5,  the  stability  index  is  plotted  as  a function  of  chord 
length  for  values  of  the  net  buoyancy  corresponding  to  16,  20,  and  24 
thousand  pounds.  The  stability  index  is  relatively  insensitive  to 
center  of  gravity  location,  hence  these  values  are  given  only  for 
3:q  = “2.0  feet.  A change  of  0.5  feet  in  the  x^  produces  no  more  than 
a 2 percent  change  in  Gn  over  the  range  of  chord  length  and  net  buoyancy 
conditions  shown,  where,  as  can  be  seen  from  Equation  (A-14),  a shift  of 
the  center  of  gravity  aft  would  tend  to  reduce  the  stability  index  and  a 
shift  forward  would  increase  the  stability  index. 

Figure  5 indicates  that  the  DOLLY  VARDEN  is  dynamically  stable  at 
infinite  speed  for  all  of  the  appendage  and  ballast  conditions  shown, 
however,  a chord  length  greater  than  approximately  5 feet  provides  little 
additional  stability. 

Also  shown  in  Figure  5 is  the  stability  index  for  the  neutrally 
buoyant  vehicle  with  the  centers  of  buoyancy  and  gravity  coincident. 

This  condition  also  provides  the  horizontal  stability  characteristics 
of  a zero  trim,  neutrally  buoyant,  60-foot  vehicle  in  powered,  horizontal 
flight. 

The  dynamic  stability  characteristics  of  the  DOLLY  VARDEN  at  finite 
speeds  will  now  be  investigated.  Figure  4 indicates  that  for  the  5-foot 
chord  length  vehicle  configuration  with  net  buoyancy  of  16,000  pounds 

and  x = -1.5  feet,  a small  positive  real  root  is  obtained  which  causes 

G 
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the  vehicle  to  be  unstable  at  finite  speeds.  Since  this  real  root 
approaches  zero  as  the  speed  becomes  large,  it  is  clear  from  Equations 
(A- 3)  through  (A-9)  and  Table  4 that  the  metacentric  effects  result 
in  the  unstable  root.  Therefore,  a variation  in  the  net  buoyancy  or  center 
of  gravity  can  alter  the  stability  characteristics  at  finite  speeds. 

The  Routh-Hurwitz  stability  criteria  (Reference  2)  for  the  cubic 
characteristic  equation  (Equation (1) ) requires  that 


A,  B,  C,  D > 0 
BC-AD  > 0 


(3) 


The  characteristic  equations  for  the  5-foot  chord  vehicle  configuration 
are  given  in  Table  4 for  the  various  ballast  conditions  indicated  in 
Table  2.  (Note  the  equations  have  been  divided  by  the  coefficient  A of 
the  cubic  term  which  is  positive  in  all  cases.)  As  can  be  seen,  for  this 
range  of  ballast  conditions  and  for  the  small  variation  in  the  hydrodynamic 
coefficients  for  the  various  chord  length  configurations  shown  in  Table  3 
the  coefficients  A,  B,  and  C are  all  positive  and  by  calculation, BC-AD>0. 
The  requirement  for  dynamic  stability  then  becomes 


D 


ZwM* 


/ I 

7‘q 


> 0 


(4) 


Substitution  of  Equations  (A-3)  and  (A-4)  into  Equation  (4)  then  provides 
the  following  condition  for  vertical  stability  at  finite  speeds  within 
the  range  of  ballast  and  appendage  configurations  considered 


L(W-B) 

XGW 


(5) 


i 
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where  L is  the  vehicle  length  and  (W-B)  is  the  weight  minus  buoyancy. 

Equation  5 indicates  that  increasing  net  buoyancy  tends  to  be  a 
destabilizing  effect  while  movement  of  the  center  of  gravity  aft 
stabilizes  the  vehicle.  The  combinations  of  net  buoyancy  and  center  of 
gravity  which  provide  dynamically  stable  motion  at  finite  speeds  are 
shown  in  Figure  6 for  three  values  of  the  appendage  chord  length.  The 
region  above  a given  curve  denotes  stability  and  region  below  denotes 
instability . 

Figure  6 indicates  that  only  marginal  improvement  in  the  range  of 
dynamically  stable  ballast  conditions  is  obtained  as  the  appendage  chord 
length  is  increased  from  5 to  7 feet.  Therefore,  based  upon  the  infinite 
speed  stability  characteristics  shown  in  Figure  5,  the  5-foot  chord 
appendage  appears  to  provide  nearly  optimum  stability  with  a minimum  of 
drag.  It  should  be  noted  from  Figure  6 that  there  will  be  certain 
ballast  conditions  for  which  the  5-foot  chord  vehicle  configuration  is 
unstable,  however,  it  will  be  shown  in  the  following  section  that  the 
instability  is  so  weak  that  the  pitch  angle  and  vehicle  trajectory  are 
not  significantly  affected  over  the  distance  of  vertical  travel. 
Furthermore,  a small  instability  would  tend  to  improve  the  response 
during  the  recovery  maneuver. 

MOTION  AND  TRAJECTORY 

In  the  previous  section,  it  was  shown  that  a 5-foot  chord  appendage 
configuration  provides  dynamic  stability  at  infinite  speed  for  the 
DOLLY  VARDEN  in  vertical  buoyant  rise  flight,  although  a very  weak 
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Pipure  6 - Effect  of  Ballast  on  Vertical  Stability 
at  Finite  Speed 
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instability  is  present  for  some  ballast  conditions  at  the  lower  speeds. 

In  this  section  the  vertical  plane  motion  of  the  5-foot  chord  vehicle 
during  vertical  flight  and  the  proposed  recovery  maneuver  is  investigated 
by  performing  computer  simulations  with  the  program  ZZMN. 

The  computer  program,  ZZMN,  is  a time  domain  simulation  essentially 
based  upon  the  nonlinear  equations  of  motion  given  in  Reference  1.  The 
coefficients  given  in  Table  3 for  the  5-foot  chord  configuration  have 
been  used  in  the  simulation  with  the  exception  that  the  static  derivatives 
Z^'  and  M^'  have  been  replaced  by  a nonlinear  fit  for  the  static  force 
and  moments,  the  coefficients  of  which  are  provided  in  Table  5.  The 
vertical  plane  equations  used  in  the  simulation  are  provided  in  Appendix  C. 

Since  the  vehicle  initially  starts  from  rest  and  accelerates  to  a 
high  speed,  the  computer  program  has  been  modified  to  include  a speed 
dependent  drag  coefficient.  The  drag  is  given  by 

Drag  = 1/2  ^(Cp  + Cr)S-j>ot^^ 


where  ST0T  is  the  total  wetted  surface  area.  The  frictional  resistance 
coefficient,  C , is  given  by  the  ITTC  curve  (Reference  3) 

r 
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F 


0.075 

(1°K,<)ReL  " 2)2 


where 


Re  - UL/w;  . 

L/ 


The  residual  resistance  coefficient,  C_ , was  estimated  for  the  5-foot 

K 

2 

chord  appendages,  and  for  STQT  = 1317  feet  , CR  = 0.000285 
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Table  5 - Nonlinear  Static  Coefficients  for 
the  DOLLY  VARDEN 

Chord  Length  » 5.0  feet 

Zw  -0.02476  ^ +0.006853 

2^  ^ -0.09560  , -0.02249 


Figures  7 and  8 show  the  velocity  and  pitch  angle,  respectively, 

attained  at  a depth  of  200  and  400  feet  as  a function  of  the  net  buoyancy. 

The  initial  depth  is  1,000  feet,  the  initial  pitch  angle  is  89  degrees,  and 

the  center  of  gravity  location  Is  1.5  feet  aft  of  the  center  of  buoyancy. 

Figure  7 shows,  as  expected,  that  an  increase  in  the  net  buoyancy  results 

in  greater  speed.  The  required  40-knot  speed  at  the  400-foot  depth  is 

obtained  with  a net  buoyancy  slightly  greater  than  20,000  pounds. 

Figure  8 shows  that  the  deviation  of  the  pitch  angle  from  the  initial 

condition  increases  with  distance  traveled  and  increasing  net  buoyancy. 

This  is  a result  of  the  weak  dynamic  instability  mentioned  in  the  previous 

section  for  the  particular  ballast  condition.  The  pitch  deviation  from 

the  vertical,  however,  is  less  than  2 degrees. 

Figure  9 shows  the  pitch  angle  and  depth  responses  during  the 

recovery  maneuver  for  a vehicle  configuration  with  net  buoyancy  of  20,000 

pounds  and  x = -1.5  feet.  The  sternplanes  are  deflected  to  15  degrees 
G 

at  a rate  of  30  degrees/second  when  a depth  of  200  feet  is  attained.  At 
the  instant  the  planes  are  deflected,  the  elapsed  time  after  launch  is 
18.9  seconds,  the  speed  is  42.4  knots,  and  the  pitch  angle  is  88.4  degrees. 
At  an  elapsed  time  of  23.2  seconds  after  launch,  the  pitch  angle  is  zero 
(body  is  horizontal).  At  this  instant,  the  depth  is  5 feet,  the  forward 
speed  is  23.4  knots,  and  the  upward  velocity  is  9.2  feet/second. 

Figure  9 shows  various  trajectory  parameters  for  the  vehicle 
calculated  with  the  initial  depth  set  in  the  simulation  at  a value  of  5000 
feet  and  the  depth  at  the  end  of  a trajectory  given  as  the  number  of  feet 
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Chord  Length  a 5.0  feet 
Center  of  Gravity  -=  -1.5  feet 
Initial  Depth  a 1000  feet 
Initial  Speed  a 0 knots 
Initial  Pitch  a 89  degrees 


Depth  a 200  feet 


Depth  > 400  feet 


16  20  24  X 103 

Net  Buoyancy  in  pounds 

Figure  7 - Effect  of  Net  Buoyancy  on  Speed  During 
Vertical  Flight 


Depth  a 400  feet 


Depth  = 200  feet 


16  20  24  X 10" 

Net  Buoyancy  in  pounds 


Figure  8 - Effect  of  Net  Buoyancy  on  Pitch  Angle 
During  Vertical  Flight 


Pitch  Angle  in  degrees 


A 


deeper  (positive)  or  more  shallow  (negative)  than  4000  feet.  The  cal- 
culated depth  change  would  be  the  same  for  any  initial  depth  selected. 

The  computer  printout  for  the  run  shown  in  Figure  9 is  provided  in 
Appendix  B. 

Figures  10  and  11  show  the  forward  speed  as  a function  of  time  for 
several  values  of  the  net  buoyancy.  In  Figure  10  the  sternplanes  are 

activated  at  a depth  of  200  feet,  in  Figure  11  at  a depth  of  400  feet. 

In  both  cases  the  vehicle  starts  from  rest  1000  feet  below  the  surface 
with  an  initial  pitch  angle  of  89  degrees.  The  sternplanes  are  deflected 
^ to  an  angle  of  15  degrees  at  a rate  of  30  degrees/second.  The  locations 

in  the  trajectory  where  the  sternplanes  are  activated  and  where  the  pitch 
angle  and/or  depth  are  zero  are  indicated. 

Since  the  vehicle  exhibits  a small  dynamic  instability  for  the 
particular  ballast  condition  considered,  the  initial  conditions  will 
i determine  how  rapidly  the  vehicle  deviates  from  a vertical  trajectory. 

, The  vertical  trajectories  of  the  DOLLY  VARDEN  for  the  above  ballast 

condition  were  also  obtained  for  two  values  of  initial  pitch  rate  with 
, an  initial  pitch  angle  of  89  degrees.  At  a depth  of  200  feet  the  pitch 

angle  was  87.2  and  85.9  degrees  for  initial  pitch  rates  of  -0.25  and 
• -0.5  degrees  per  second,  respectively.  Computer  printouts  of  these 

runs  are  also  provided  in  Appendix  B. 

The  effect  of  maximum  sternplane  angle  and  rate  on  the  recovery 
maneuver  has  also  been  investigated.  Key  parameters  of  the  trajectory 
are  shown  in  Tables  6 and  7,  respectively,  for  maximum  sternplane  angles 
v from  12  to  21  degrees  deflected  at  a rate  of  30  degrees/second,  and  stern- 

plane rate  of  7.5  to  75  degrees/second  deflected  to  a maximum  angle  of 
15  degrees.  The  planes  are  deflected  at  a depth  of  4200  feet  or  18.9 
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Sternplane  Deflection  to 
15  degrees  at 
30  degrees/second 

^Initial  Pitch  « 89  degrees 


Figure  10  - Velocity  Trajectories  for  Recovery 
Maneuver  Initiated  at  200  Feet 
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Speed  in  knots 


Sternplane  Deflection  to 
15  degrees  at 
30  degrees/second 
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Table  6 - Trajectory  Parameters  for  Sternplane 
Deflection  Angle  Variation 

Deflection  Rate  - 30  Degrees/second 


Tra  je 

ctor  y 

Pa  r a meters 

Elasped  time  in  seconds 

f rom 

launch 

to : 

Ze  ro 

Pitch 

Ze  ro 

De  pt  h 

Speed 

in  knots  at 

Ze  ro 

Pitch 

Ze  ro 

Dept  h 

Do  pt  h 

in  fe 

et  at  : 

Ze  ro 

Pitch 

Pitch 

in  degrees 

at  : 

Zero 

Dept  h 

Upwa  r 

d vert 

ica  l 

vo  1 or 

i t y in 

feet 

sec  at  : 

Ze  ro 

Pitch 

Hor  iz 

onta  I 

distance 

t rave 

led  in 

feet 

a t : 

Zero 

Pi  t ch 

Zero 

De  pt  h 

Maximum 

Sternplane  Deflection  Angle 

in  degrees 

12 

15 

18 

21 

23  .7 

23 .2 

22.9 

22  .7 

22.4 

23.4 

““ 

25.6 

23  .4 

21  .4 

19.7 

30.6 

22.8 

-22.0 

1 .4 

17.0 

28.0 

20.1 

-2.5 

- 

8.8 

9.2 

9.5 

9.7 

147 . 5 

123.8 

107.2 

94 .9 

92.9 

129.7 
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Table  7 - Trajectory  Parameters  for  Sternplane 
Deflection  Rate  Variation 

Deflection  Angle  * 15  Degrees 


Trajectory 

Sternplane  Deflection  Rate 
in  degrees/second 

ra  r a in©  ters 

m 

15 

30 

50 

75 

Elasped  time  in  seconds 
from  launch  to: 

Zero  Pitch 

23.5 

23  .4 

23.2 

23 .2 

23 .1 

Zero  Depth 

21  .9 

22  .4 

23.4 

Speed  in  knots  at: 
Zero  Pitch 

24  .6 

23.7 

23  .4 

23  .3 

23  .3 

Zero  Depth 

34 .9 

28.2 

22  .8 

— 

Depth  in  feet  at : 
Zero  Pitch 

-47  .6 

-15.5 

1 .4 

8.3 

11.8 

Pitch  in  degrees  at: 
Zero  Depth 

38.6 

16.3 

-2.5 

- 

- 

Upward  vertical 
velocity  in  feet/sec  at: 
Zero  Pitch 

9.4 

9.3 

9.2 

9 .2 

9.2 

Horizontal  distance 
traveled  in  feet  at: 
Zero  Pitch 

129.0 

124.8 

123.8 

123  .6 

123 .5 

Zero  Depth 

53 .1 

86.7 

129.7 

30 


42.4 


seconds  after  launch,  at  whic  ^ime  the  initial  conditions  are  “ 
knots  and  the  pitch  is  88.4  degrees.  A blank  space  indicates  the  vehicle 
did  not  reach  a depth  of  4000  feet  after  being  released  from  5000  feet. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A cruciform  configuration  of  all-movable,  rectangular  stern  append- 
ages having  a tip  to  tip  span  of  10  feet  and  chord  length  of  5 feet 
provides  adequate  stability  and  maneuverability  for  operational  require- 
ments of  the  60-foot  DOLLY  VARDEN. 

1.  The  vehicle  exhibits  a very  weak  pitch  instability  for  some 
ballast  conditions  in  vertical,  buoyant  rise  flight,  however,  the 
instability  is  so  weak  that  the  pitch  deviation  from  the  vertical  is 
less  than  5 degrees  for  an  initial  pitch  angle  of  89  degrees  and  initial 
pitch  rate  of  5 degrees  per  second. 

2.  With  a net  buoyancy  of  20,000  pounds,  and  the  center  of  gravity 
location  1.5  feet  aft  of  the  center  of  buoyancy,  the  required  speed  of 
40  knots  is  attained  at  the  400-foot  depth. 

3.  Recovery  from  vertical  flight  at  a depth  of  200  feet  can  be 
achieved  by  deflection  of  the  proposed  sternplanes  to  15  degrees  at  a 

rate  of  30  degrees/second.  This  maneuver  brings  the  vehicle  to  a horizontal 
attitude  near  the  surface  with  a forward  speed  of  approximately  23  knots  and 
an  upward  vertical  velocity  of  9.2  feet/second. 

ACKNOWLEDGMENTS 

The  author  wishes  to  acknowledge  the  contributions  of  Elizabeth 
Dempsey  and  William  Day  of  the  Ship  Performance  Department  for  their 
estimation  of  the  hydrodynamic  and  drag  coefficients  for  the  vehicle. 

He  also  wishes  to  thank  Dr.  J.P.  Feldman  for  his  valuable  suggestions 
and  review  of  the  material  presented  herein. 


REFERENCES 


- 


1.  Gertler,  M.  and  G.R.  Hagen,  "Standard  Equations  of  Motion  for 
Submarine  Simulation,"  DTNSRDC  Report  2510,  June  1967. 

2.  Jones,  A.L.  and  B.R.  Briggs,  "A  Survey  of  Stability  Analysis 
Techniques  for  Automatically  Controlled  Aircraft,"  NACA  Technical 
Note  2275,  January  1951. 

3.  Principals  of  Naval  Architecture,  Society  of  Naval  Architects  and 
Marine  Engineers,  1967,  p.  298. 


J 


i- 


32 


! 


APPENDIX  A 

EQUATIONS  FOR  DYNAMIC  STABILITY' 


Y\ 


1 


33 


APPENDIX  A 


EQUATIONS  FOR  DYNAMIC  STABILITY 

The  equations  for  the  dynamic  stability  of  the  DOLLY  VARDEN  buoyant 
rise  vehicle  are  derived  in  this  section.  Linear  equations  of  motion 
in  the  vertical  plane  are  formulated  for  small  perturbations  about  an 
arbitrary  equilibrium  body  attitude  and  speed,  and  the  net  buoyancy  and 
longitudinal  center  of  gravity  location  are  retained  as  arbitrary 
variables.  Because  of  the  geometric  symmetry  of  the  DOLLY  VARDEN  relative 
to  the  vertical  and  horizontal  planes,  the  vertical  plane  equations  are 
also  applicable  to  the  horizontal  plane  when  the  net  buoyancy  is  zero 
and  the  centers  of  buoyancy  and  gravity  are  coincident. 

The  formulation,  notations,  and  conventions  of  this  section  follow 
those  commonly  used  in  the  standard  equations  for  submarine  motion  as 
given  in  Reference  1.  The  equations  of  motion  are  formulated  relative 
to  a body  fixed  Cartesian  coordinate  system  with  the  origin  located  at 
the  center  of  buoyancy.  The  positive  x axis  is  directed  forward  along 
the  axis  of  symmetry,  and  the  positive  z axis  is  directed  downward  when 
the  vehicle  is  in  a horizontal  attitude. 

In  the  equilibrium  condition  the  vehicle  is  assumed  to  be  moving  in  a 

purely  forward  direction  with  constant  speed  U at  an  angle  6o  relative  to 

the  horizontal  (0  = n/2  corresponds  to  vertical  flight.)  It  should  be 

o 

emphasized  that  the  DOLLY  VARDEN  buoyant  rise  vehicle  will  rapidly  accel- 
erate upward  initially,  thus  not  strictly  satisfying  the  equilibrium 
steady  speed  condition.  Therefore,  the  stability  analysis  described  herein 
better  represents  the  higher  speed  than  the  lower  speed  flight  of  the  vehicle. 
The  body  is  perturbed  in  the  pitch  mode  and  the  subsequent  vertical  plane 
motions  (x-z  plane)  are  to  be  determined.  The  perturbation  velocities  in  the 


x and  z directions  and  the  pitch  perturbation  angular  velocity  are 
denoted  by  u,  w,  and  q,  respectively.  Due  to  gravitational  effects, 


displacement  dependent  terms  appear  in  the  equations  of  motion,  and 


the  perturbation  pitch  displacement  is  defined  by 


*=/?«*  • 


Under  the  condition  that  the  center  of  gravity  of  the  body  is 
located  on  the  body  axis  of  symmetry,  the  equations  for  the  normal  force 
and  pitch  moment  may  be  decoupled  from  the  axial  motion  to  first  order 
in  the  perturbation  quantities.  The  normal  force  and  pitch  moment  are, 
respective ly . 


- f/zlH  Zfg  + %.LJZ*w 

•f  fa  L U +■  i fa  L V 0 (a-i; 


Ljg  = fa  w 

+ faC  v /n'gf  + fa l ti/n'wv  + fai? 


(A-2 ) 


where  the  hydrodynamic  coefficients  are  given  in  nondimensional  form  as 

defined  in  Reference  1.  The  metacentric  derivatives,  Z ' and  M ' , are 

u u 

given  about  the  equilibrium  pitch  angle  0o  as  follows: 


z*  = 
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For  this  investigation  the  values  of  z^,  z^,  and  Xg  are  zero.  In  addition, 


Hence,  Equations  (A-l)  and  (A-2)  form  a pair  of  linear,  homogeneous, 
differential  equations  for  w and  6.  The  nondimensional  roots  o^'  of 
the  characteristic  equation  for  this  set 

+ Bo;1  -I-  C<r'  + D = 0 

provide  the  basic  form  of  the  solution  which  is 


K-l 


where : 
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and 
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A stable  solution  for  0 and  w Is  obtained  when  all  of  the 
roots  of  Equation  (A-5)  have  negative  real  parts.  The  Routh-Hurwitz 
criterion,  as  outlined  in  Reference  2,  provides  a simple  test  for 
stability.  A necessary  and  sufficient  condition  for  all  the  roots  of 
Equation  (A-5)  to  have  negative  real  parts  is  that 


A,  B,  C,  D > O 
pc  -Ad  >0 


(A-10) 


At  infinite  speed,  the  nondimensional  metacentric  derivatives 
and  Me  are  zero  and  the  last  term,  D,  in  Equation (A-5)  is  zero.  One 
root  of  Equation  (A-5)  is  then  zero  and  the  other  two  roots  are  obtained 
by  the  solution  of  the  quadratic. 


I 


A<rz  + Bcr  + C0  = o 


(A-ll) 


where  A and  B are  given  from  Equations  (A-6) and  (A-7),  respectively, 

and  C is  defined  by, 
o 

c = zl  ( m'j  -/»'<- ) - /n'u  ( ~£t  + ni ) <*-“) 

A necessary  and  sufficient  condition  that  the  roots  of  Equation  (A-ll) 
have  negative  real  parts  is  that 
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.1 


For  most  body  shapes  having  a reasonable  degree  of  longitudinal 
symmetry  and  uniform  mass  distribution,  the  coefficients  A and  B given 
by  Equations  (A-6)  and  (A-7)  are  positive.  The  infinite  speed  stability 
criterion  Equation  (A-13)  then  reduces  to  C0  ^ O .A  measure  of  the 


stability  of  the  system  can  be  given  in  terms  of  the  stability  index  G^, 


where  for  stability 


n r,  - i - ■ ■*- — - — 


and  an  increasing  value  of  G^  indicates  a greater  degree  of  stability. 


HORIZONTAL  FLIGHT 

With  the  exception  that  gravitational  terms  are  not  present  in  the 
horizontal  plane,  the  linearized  stability  equations  for  coupled  yaw  and 
lateral  motion  are  analogous  to  the  vertical  plane  equations  at  infinite 
speed  where  in  Equations  (A-l)  and  (A-2)  v and  r,  the  perturbation  veloc- 
ities in  the  lateral  and  yaw  modes, replace  w and  q,  respectively,  and  the 
horizontal  plane  derivatives  Y’  and  N1  replace  the  corresponding  vertical 
plane  derivatives  7.'  and  M'  , respectively,  as  indicated  in  Table  4. 

The  DOLLY  VARDEN  is  geometrically  symmetric  relative  to  the  hori- 
zontal and  vertical  planes,  thus  the  characteristic  equation  and  stability 
criteria  corresponding  to  infinite  speed,  vertical  plane  motion  (Equations 
(A-ll)  through  (A-14))also  provide  the  stability  conditions  for  the  hori- 
zontal plane. 
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EQUATIONS  USED  IN  TRAJECTORY  SIMULATION 


The  non-linear  equations  used  to  simulate  the  vertical  plane  trajec- 
tory of  the  DOLLY  VARDEN  are  the  "Standard  Equations  of  Motion  for  Submarine 
Simulation"  given  in  Reference  1.  In  the  present  application  all  terms 
involving  horizontal  plane  variables  (v,  r,  and  p)  are  zero  and  the  roll 
angle  is  zero.  In  addition,  some  non-linear  vertical  plane  coefficients 
were  not  estimated,  and  their  effect  in  the  equations  is  neglected.  The 
axial  force,  normal  force,  and  pitch  moment  equation  are  given  below 
where  the  hydrodynamic  coefficients  are  given  in  Tables  3 and  5 and  X ' 
is  given  bv  the  equation 


* ~ ( C-r  + ) U* 


Axial  Force 

1 + = 

f/z  L*  'xl*  uf  - ( w - B ) s \n.e 


m (u  + ufg  ) - m y 


(C-l) 


(C-2) 


Normal  Force 


m(ih  - ng  ) - = (c- 

+ ^ + ^ ^ ^ ^ * fa L ^ 
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Pitch  Moment 


Ijj  ^ + ”i-j6  (-9.  t - tnxc.  ( - ufl  ) — 

+ h ls  m'j  i + fo  l'  /n*  w + f/j.  l*  /r£  «? 

+ ^ L3  [ <,  UW  -f  /MW/UJ/  TA/lu\~\ 

-(X^W-XgB)  &>s$  “ ( Jc-W B) 


(C-4) 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

(1)  DTNSRDC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

12)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS.  NOT  FOR  GENERAL  DISTRIBUTION 


